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ABSTRACT 

The  Mode-I  fracture  toughness  of  a  brittle  material  reinforced  by  aligned  brittle  fibers  is  studied 
theoretically.  The  fibers  are  assumed  to  slip  relative  to  the  matrix  when  a  critical  interface  shear  stress  is  reac'.ed, 
and  the  toughening  action  of  the  fibers  is  presumed  to  be  due  to  bridging  of  crack  faces  in  the  vicinity  of  the  crack 
front.  The  toughening  due  to  the  fiber  reinforcement  is  related  to  basic  parameters  associated  with  the  related 
problem  of  steady-state  matrix  cracking  in  the  presence  of  intact  fibers.  Bridge  lengths  at  fracture  and  fracture 
resistance  curves  are  calculated.  \  '  •  •  '  \  ^  M  \  ■  \  ;'  r  •  t  ^  o\ ,  *  -  '  L-  » 

NOMENCLATURE 


a  fiber  radius 

c  fiber  volume  fraction 

E  longitudinal  composite  Young's  modulus,  =  cEf -t- (l-c)Eni 
Ef,  Em  fiber,  matrix  Young's  moduli 
Gf,  Gm  fiber,  matrix  shear  moduli 

Gm  critical  matrix  energy-release-rate,  =  (l-v^mJK^m/Em 
K  applied  stress-intensity  factor 
Km  critical  matrix  stress-intensity  factor 

L  bridge  length 

Q  complementary  energy,  spring  model 
S  fiber  breaking  stress 

V  strain  energy,  spring  model 

V  crack-face  displacement 

a  non-dimensional  bridge  length,  Eq.  (26) 

A  modified  toughening  ratio,  Eq.  (19) 


*  On  leave  from  the  Department  of  Mathematics,  University  of  Nigeria,  Nsukka,  Nigeria. 


X  toughening  ratio,  K/Km 
Vf,  Vm  fiber,  matrix  Poisson's  ratio 
ao,  Cl  reference  stresses,  Eqs.  (1-2) 

Ccr  critical  stress,  steady-state  matrix  cracking 
T  fiber-matrix  interface  shear  resistance 


INTRODUCTION 

'^his  paper  addresses  the  problem  of  calculating  the  increase  in  the  Mode  I,  plane-strain 
fracture  toughness  of  a  brittle  material  when  it  is  reinforced  by  long,  aligned,  brittle  fibers,  with 
sliding  between  the  fibers  and  the  matrix  suppressed  only  if  the  interface  frictional  shear  is  less 
than  some  limiting  stress.  The  configuration  contemplated  (Fig.  1)  is  an  infinite  domain  containing 
a  semi-infinite  crack  that  is  bridged  by  intact  fibers  in  the  vicinity  of  the  crack  tip.  The  crack  is 
growing  in  a  quasi-static,  steady-state  fashion,  with  simultaneous  fracture  of  the  matrix  along  the 
crack  front  and  failure  of  the  fibers  at  the  end  of  the  bridged  zone.  This  crack  propagation  is 
imagined  to  occur  under  the  imposition  of  a  remote  stress  field  that  corresponds  to  an  "applied" 
stress  concentration  factor  K,  which  is  accordingly  defined  as  the  fracture  toughness  of  the 
composite  material.  The  primary  aim  of  the  present  study  is  to  provide  theoretical  results  for  the 
:  Highening  ratio  X=K/K,n,  where  K„  is  the  fracture  toughness  of  the  unreinforced  matrix  material. 
In  addition,  resistance  curves  will  be  produced,  showing  how  an  initially  unbridged  crack  grows 
into  the  matrix  material  as  the  applied  K  is  increased  towards  its  critical  value. 

The  present  results  should  be  applicable  to  less  idealized  geometries  if  the  conditions  of 
small-scale  bridging  are  met,  wherein  the  bridge  length  L  is  small  relative  to  all  other  pertinent 
dimensions,  such  as  crack  length  and  distance  of  the  crack  tip  from  the  boundary.  Work  on  the 
toughening  problem  by  MarsH/XL  and  EVANS  (1986)  and  MARSHALL  and  COX  (1987) 
considered  configurations  for  which  the  small-scale  assumption  was  not  made.  Besides  its 
emphasis  on  small-scale  bridging,  the  present  paper  differs  from  these  studies  primarily  in  the  way 
the  analysis  and  results  are  linked  closely  to  basic  parameters  associated  with  steady-state  matrix 
cracking  without  fiber  failure  (BUDIANSKY,  HUTCHINSON,  and  EVANS,  1986).  The  importance 
of  a  statistical  variation  in  fiber  strength  have  recently  been  emphasized  by  THOULESS  and  Evans 
(1988),  and  the  influence  of  residual  stresses  has  been  discussed  by  MARSHALL  and  Evans 
(1988),  but  these  effects  are  not  considered  in  the  present  paper,  which  is  based  upon  the 
preliminary  analysis  that  was  outlined  by  BUDIANSKY  (1986). 
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Steady-State  matrix  Cracking 


We  review  here  results  found  by  BUDIANSKY,  HUTCHINSON,  and  EVANS  (BHE)  for  the 
stress  Ocr  associated  with  steady  matrix  cracking  (Fig.  2).  In  the  absence  of  initial  stresses,  Ocr  is 
given  in  Fig.  3  in  terms  of  the  reference  stresses  ctq  and  a\  defined  by 
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where 
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The  ratio  Oi/oq  plays  a  key  role  as  a  non-dimensional  parameter  that  characterizes  the  frictional 
shear  strength  of  the  fiber-matrix  interface.  For  x  0,  Oi/oq  0>  and  then  Ocr  Ofi  thus,  Oi  is 
the  critical  stress  for  matrix  cracking  when  there  is  extensive  interface  slip  between  the  matrix  and 
the  bridging  fibers.  On  the  other  hand,  for  Oi/oq  >  there  is  no  interface  slip  during  steady- 
state  matrix  cracking,  and  in  this  case  Ocr  =  Oq. 

For  Oi/oq  <  3^^3^  the  curve  in  Fig.  3  was  described  parametrically  in  BHE  by 
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for  1  <  Y  <  oo.  This  provides  the  implicit  equation 
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for  Ocr- 

If  the  fibers  remain  intact  after  the  matrix  is  fully  cracked,  the  composite  may  continue  to  , 
resist  additional  stress,  as  shown  schematically  by  the  idealized  stress-strain  relation  in  Fig.  4.  The 
ultimate  composite  strength  is  then  approximately  Omax=  cS,  where  S  is  the  breaking  stress  of  the 
fibers.  .;.i 
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Spring  Model;  Equiyalent  Nonlinear  Spring  Stiffness 


The  steady-state  matrix  cracking  problem  of  Fig.  2  may  be  modeled  by  the  configuration 
shown  in  Fig.  5a,  wherein  crack-bridging  springs  of  zero  initial  length  are  imagined  to  emerge 
continuously  from  the  crack  tip  as  the  crack  grows.  The  stress  o  in  the  springs  will  be  assumed  to 
depend  on  the  the  crack  opening  5=2v  (Fig.  6)  in  such  a  way  as  to  make  the  matrix  cracking  stress 
of  the  model  the  same  as  that  of  the  fibrous  composite.  In  establishing  this  equivalence,  the  crack- 
tip  energy-release-rate  in  the  model  problem  will  be  modified  by  the  factor  (1-c)  in  order  to 
simulate  the  reduction  in  crack  area  associated  with  the  presence  of  fibers.  Once  the  appropriate 
spring  characteristics  are  thereby  deduced  from  the  SHE  results  for  steady-state  matrix  cracking, 
we  will  be  able  to  use  the  spring  model  to  study  the  toughening  problem  (Fig.  1). 

By  setting  the  modified  crack-tip  energy-release-rate  equal  to  the  difference  between  the  far 
upstream  and  downstream  potential  energies  per  unit  length,  we  get  the  condition 


od5  =  (l-c)G„, 


(6) 


governing  the  critical  matrix-cracking  stress  of  the  model  problem.  Alternatively,  we  can  write  the 
J-integral  (RICE.  1^68) 


J  = 


dup 
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around  the  path  shown  in  Fig.  5b  to  get  the  same  result,  if,  again,  we  multiply  the  contribution 
from  the  small  circle  around  the  crack  tip  by  the  factor  (1-c).  Eq.  (6)  is  equivalent  to 


Q(Ocr) 


where  (see  Fig.  6)  Q  is  the  complementary  energy  function 


(8) 
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Note,  for  future  reference,  that  the  strain  energy 


is  related  to  Q(a)  by 


V(c)  =  aQ’(a)  -  Q(a). 


We  now  demand  that  the  cracking  criterion  (8)  give  the  BHE  results  for  Ocr-  In  the  slipping-fiber 
case,  we  equate  the  left-hand  sides  of  (5)  and  (8)  to  get 


Q(CTcr) 


laiJ  SvOiJvaoJ  H\Cq) 


In  the  no-slip  case,  the  cracking  criterion  may  be  written  as  (Ocr/Oo)^  =  1.  and  comparison  of  this 
with  (8)  gives 


Q(CJcr) 

(l-c)G„ 


(The  validity  of  these  identifications  is  corroborated  by  the  fact  that  the  right-hand  sides  of  both 
(12)  and  (13)  are  indeed  proportional  to  l/Gm-)  We  can  now  write  a  general  formula  for  Q(o)  by 
noting  that  the  no-slip  condition  for  fibers  (Eq.  (27)  of  BHE)  may  be  expressed  as 
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in  terms  of  the  smeared-out  fiber  stresses  0  at  the  crack  face.  Accordingly,  the  complementary 
energy  function  for  the  springs  is 


^0 


From  (11),  the  strain  energy  function  is 
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(This  last  equation  for  V  corrects  errors  in  Eq.  (28)  of  BUDIANSKY  (1986)).  Finally,  the  stress- 
displacement  relation  for  the  springs  is  given  by  v  =  5/2  =  Q'(a)/2,  whence 
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Fracture  toughness 

We  now  contemplate  the  configuration  in  Fig.  1,  and  model  it  with  that  shown  in  Fig.  7a, 
wherein  the  semi-infinite  crack  is  partially  bridged  by  springs  possessing  the  properties  that  we 
have  just  established.  Far  from  the  crack  tip  the  material  is  regarded  as  homogeneous  but 
orthotropic,  and  subjected  to  the  far-field  stress  state 


V27tr 


fap(e) 


where  fap(O)  =  1,  and  the  angular  distribution  function  fa3(6)  =  fap(-6)  is  appropriate  for  the 
particular  orthotropy  of  the  composite.  Thus  K  is  an  "applied"  Mode-I  stress-intensity  factor, 
while  the  stress-intensity  factor  along  the  crack  edge  in  the  matrix  material  is  assumed  equal  to  Km- 
Now  write  the  modified  J-integral  around  the  path  in  Fig.  7b  to  get 
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+  V[a(L)] 


where  the  constant  A  takes  into  account  the  orthotropy  of  the  composite,  and,  as  before,  the  (1-c) 
factor  reflects  the  fact  that  fibers  interrupt  the  crack  front  Note  that  the  first  term  on  the  right  is  the 
same  as  (l-c)Gm.  If  we  now  assert  that  the  bridged  configuration  propagates  steadily,  with 
simultaneous  matrix  fracture  and  "last"  fiber  failure,  we  can  set  a(L)  equal  to  cS,  where  S  is  the 
fiber  strength.  Now  define  the  modified  toughening  ratio 

K/K„ 

V(l~c)AE/Em 

and  use  the  formulas  (16)  for  V  in  (18).  Then  we  get 
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and  this  gives  the  curves  in  Fig.  8,  showing  A  vs.  cS/Oi  for  Oi/oq  =  0,1,2.  An  estimate  for  the 
orthotropy  factor  A  in  the  definition  (19)  of  the  modified  toughening  ratio  A  is  given  in  Appendix 
A. 

These  results  may  be  rewritten  to  show  A  as  a  function  of  (cS/Ocr)  in  the  form 
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wherein  (Ocr/Oo),  given  in  Fig.  3,  in  turn  depends  on  Oi/oq.  Thus,  the  toughening  has  been 
related  directly  to  the  ratio  of  the  fiber  strength  to  the  theoretical  stress  for  steady-state  matrix 
cracking.  Curves  showing  A  vs.  cS/Ocrfor  (ai/ao)^  =  0,  1,  2,  3,  and  «»  are  given  in  Fig.  9. 
Recall  that  in  the  case  of  steady-state  matrix  cracking,  the  post-cracking  strength  amax  of  the 
composite  is  approximately  equal  to  cS,  assuming  cS>Ocr.  Accordingly,  the  abscissa  in  Fig.  9 
may  be  interpreted  as  Omax/Ocr  for  cS/Ocr  >1. 

In  the  limiting  case  Oi/ag  — > «»,  Ocr  =  Oq,  and  the  modified  toughening  ratio  is  simply 
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(For  all  finite  Oi/ao  >  3^^^^  vve  still  have  Ocr  =  Oq,  but  then  (22)  holds  only  for  cS/oq  <(ai/ao)^/3; 

for  larger  values  of  cS/Oq  frictional  sliding  will  occur  before  fracture,  even  if  it  is  absent  during 
steady  matrix  cracking.) 

For  oi/ao  — >  0,  in  which  case  (Jct  =  Oi,  we  have 
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Fig.  8  shows  that  this  approximation  remains  quite  accurate  for  moderately  high  values  of  ai  /OQ. 

Note  that  Figs.  3  and  9  show  that  for  decreasing  frictional  resistance  T,  and  hence  lower  Oi, 
the  steady-state  matrix-cracking  stress  Ccr  Soes  down,  but  the  fracture  toughness  increases.  The 
curves  in  Fig.  9  would  appear  to  imply  an  inherent  limitation  on  the  amount  of  toughening  that 
could  reasonably  be  expected  from  a  well-designed  composite.  Fracture  toughness  is  an  increasing 
function  of  the  ratio  associated  with  matrix  cracking  and  failure  in  uniaxial  tension.  Since 

it  is  unlikely  that  values  much  larger  than  2  for  this  ratio  would  be  considered  desirable, 
toughening  ratios  around  4  may  be  the  most  one  might  seek. 

TOUGHENING  VS.  BRIDGE  LENGTH;  RESISTANCE  CURVES 

To  calculate  the  size  of  the  bridged  zone  when  fracture  occurs  we  have  to  set  up  and  solve 
an  integral  equation  for  the  stress  in  the  intact  fibers  along  the  crack  faces.  For  this  purpose  we 
contemplate  the  spring  model  of  Fig.  7a,  and  note  that  the  crack  opening  displacement  v(x)  is  given 
by 


v(x)  = 
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where  the  first  term  represents  the  displacement  due  to  the  far-field  loading,  and  the  second  term  is 
the  crack-closing  effect  of  the  spring  stresses.  The  orthotropy  factor  A,  the  same  one  that  appears 
in  the  J-integral  (18),  correctly  modifies  standard  formulas  for  crack-face  displacement  in  isotropic 
materials  (e.g.  Tada  et  al.,  1985).  Equating  the  displacements  in  the  bridged  zone  to  those  given 
by  the  constitutive  relations  (17)  found  for  the  springs  provides  the  integral  equation 
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It  is  shown  in  Appendix  B  how  this  integral  equation,  suitably  non-dimensionalized,  was  solved 
simultaneously  with  the  J-integral  relation  (18)  to  provide  results  for  the  modified  toughening  ratio 
A  as  a  function  of  the  non-dimensional  bridge  length 


=  l_E2_ffS£  1 T 

7CAE(l-C)UmJ  (26) 

Fig.  10  shows  curves  thereby  calculated  for  A  vs.  a,  for  (oi/oo)^  =  0,  2, «».  Remarkably,  the 
introduction  of  the  theoretical  steady-state  matrix-cracking  stress  Ocr  into  the  definition  (26)  for  the 
non-dimensional  bridge  length  has  made  the  results  come  out  nearly  the  same  over  the  full  range  of 
composite  parameters. 

With  A  given  by  Eq.  (21)  (or  Fig.  9),  these  results  provide  the  bridge  length  L 
corresponding  to  propagation  of  the  cracked  configuration  of  Fig.  1,  with  simultaneous  matrix 
fracture  and  fiber  failure  at  the  end  of  the  bridged  zone.  However,  the  curves  of  Fig.  10  have 
another  interpretation  as  resistance  curves  ,  in  the  following  sense.  Suppose  a  pre-existent, 
unbridged  crack  in  the  composite,  cutting  through  both  fibers  and  matrix,  is  subjected  to  a 
gradually  increasing  far-field  K,  and  re-define  A  in  Fig.  10  in  terms  of  this  current  value  of  K.  As 
K  increases,  the  crack  will  advance  into  the  matrix  material,  and  crack-bridging  fibers  will  remain 
intact  as  long  as  A  remains  below  its  critical  value  for  final  fracture  of  the  composite.  If  we  now 
redefine  L  as  the  amount  of  matrix  crack  growth  that  occurs  before  overall  composite  fracture,  the 
relation  between  the  current  values  of  K  and  L  is  given  by  the  curves  of  Fig.  10. 

It  is  noteworthy  that  reduction  of  the  frictional  resistance  x  ,  and  hence  Ocr,  will  result  in 
increased  matrix  cracking  and  bridge  lengths  before  and  at  final  fracture . 

Concluding  Remarks 

The  fracture  toughness  of  brittle  materials  reinforced  by  aligned  brittle  fibers  has  been 
related  simply  to  parameters  associated  with  the  steady-state  matrix  cracking  of  such  materials.  In 
each  of  these  processes  the  fiber-matrix  interface  shear  resistance  plays  an  important  role,  with 
opposite  effects:  decreasing  shear  strength  increases  the  toughness  but  lowers  the  matrix  cracking 
resistance.  This  is  consistent  with  the  fact  that  low  interface  friction  may  result  in  large  amounts  of 
matrix  cracking  associated  with  the  fracture  resistance  curves  of  the  composite.  The  effects  of 
various  other  physical  and  geometrical  parameters  may  be  assessed  with  the  help  of  the  formulas 
and  curves  that  have  been  presented. 
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APPENDIX  A 


Orthotropic  Factor  a 


Consider  an  orthotropic  material,  isotropic  in  planes  normal  to  the  xi-axis.  The  energy 
release  rate  Gj  associated  with  a  mode  I  crack  in  the  X2-X3  plane  is  related  to  the  stress  intensity 
factor  Kj  by  (Tada  et  al  1985) 

Gj  =  C  K?  (Al) 

where  C  is  defmed  in  terms  of  compliances  Ajj  by 


A  _  7^11^22  7^22 

^  ■  V  2  [V  All 
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2Ai2  +  A66 


-1 1/2 
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(A2) 


The  compliances  may  be  expressed  in  terms  of  conventional  elastic  constants  E,  E,  G,  v  and  v 
defined  by 

ei=Oi/E  ,  E2  =  -vai/E  for  02  =  03  =  0  ; 

£2  =  O2/E  ,  £3  =  -  VO2/E  for  o,  =  O3  =  0  ; 


and 

For  plane-strain  cracks 


Yi2  =  0.2/G. 


1 


1 


All  -  gCl  “  V^E/E),  Ai2  -  -  g  (1  v),  A22  ~  g  A66  =  Q 

In  th**  J-integral  expression  (18)  C  is  represented  by  (l-v^)/AEn,and  therefore 

1-vi 


(A3) 


A  =  ^ 


CE 


(A4) 
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For  aligned  fibers  the  v;lastic  constants  E,  E,  G,  v  and  v  may  be  estimated  on  the  basis  of  the  self- 
consistent  model  (Hill,  1965)  in  terms  of  Em,  Ef,  Vm,  Vf  and  c.  For  Vf  =  Vm  =  1/4,  the 
consequent  dependence  of  A  on  c  is  plotted  for  various  values  of  Ef/Em  >  1  in  Figure  1 1 . 


APPENDIX  B 


INTEGRAL  EQUATION  SOLUTION 


Tlie  substitutions 
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into  the  integral  equation  (25)  give 
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VT  -  f  g(t)  log  f  I  dt  =4^  g(s), 
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where 

r  =  CTi/ao,  4  =  CTc/cTo. 
Making  the  same  substitutions  (Bl)  into  (18)  gives 
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An  alternative  expression  for  A  that  can  be  shown  to  be  consistent  with  (18)  and  the  integral 
equation  (25)  follows  from  the  relation  (Tada  et  al,  1985) 


namely, 


A  = 


Jo  V7 


ds 
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(B4) 


The  solution  of  the  integral  equation  (B2)  follows  the  set-up  developed  in  BUDIANSKY, 
AMAZIGO,  and  Evans  (1985?)  Differentiate  the  equation  with  respect  to  s  to  get 
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where  Sq  satisfies 


8(so)  = 


To  satisfy  the  original  equation  (B2)  we  have  to  enforce  it  at  one  point  in  addition  to  satisfying 
(B5)  for  s  in  (0,a).  Enforcing  the  integral  equation  at  s  =  a  gives 
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We  now  let 


in  (B5)  to  obtain 
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g(s)  = 


f  —  dt  =  ^t^VT  [2sH’(s)  -  H(s)] ,  s<So 

J  0  ?  “  * 

=  ^^H^H(s)[2sH’(s)-H(s)],  s>So. 


It  is  convenient  to  write 


H(s)  =  C(s/a)  +  f(s) 


(B9a) 
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where  C  =  V^g(a),  and  f(0)  =  f(a)  =  0.  Then 
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3V2a(^]  [(^]c{C+2af(s)}  +  f{2sf'(s)-f(s)}  ]  ,  s>So  (Bllb) 


Now  let 


and  expand  f  in  the  Fourier  series 


S  =  (1  -  cos  0) 


f(s)  =  an  sin  (nG). 
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(B13) 


We  note  that 


a 

Jo  ^  X  ^ 

We  now  substitute  the  above  into  (B1 1),  multiply  the  resulting  equation  by  sin(m0)  and  integrate 
over  [0,a]  to  get 

N  N 

3  3 

B^CGo)  C  +  3VTa(^)  DJ0o)  +  X  Etnn(0o)  an  +  3V2 a(^)  cXnFnm(eo)  an 

n=l  n=l 

N  N 

+  3VTa(^)  XX^ninj(eo)  anaj  =  An,,  m=l,...N  (B14) 


where 
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■,  for  n  9^  j+m  or  j  n+m  or  m  ^  j+n 


m  =  j  +  n 


and  similar  expressions  for  n  =  j  +  m  and  j  =  n  +  m  obtained  by  cyclic 
permutation  of  j,  n  and  m. 
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Lnjn  =  sin  (0/2)  tan  (0/2)  cos  (n0)  sin  (m0)  d0, 
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were  evaluated  numerically  in  spite  of  the  fact  that  the  last  two  integrals  can  be  evaluated 
analytically  in  terms  of  finite  series.  0o  satisfies  the  equation  (B6)  with  so  =  (a/2)(l  -  cos  0o),  that 


1  r^'/a 

-(l-cos0o)C  +  2^  aBsin(n0o)  -  sin  (0o/2)  =  0. 


(B15) 


Substitution  for  g  into  (B7)  using  (B8),  (BIO),  (B12)  and  (B13)  gives 
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Equation  (B3)  relating  C  and  A  is 
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while  the  alternative  formula  (B4)  becomes 


A  = 
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-  C  - 


an 
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(B18) 


For  prescribed  values  of  r  =  Oi/Ooand  a  the  N+3  nonlinear  equations  (B14),  (B15),  (B16)  and 
(B17)  were  solved  by  the  Newton-Raphson  iteration  method  for  the  N+3  unknowns  C,  an, 
n=l,..,N,  A  and  Gq.  The  formula  (B18)  provided  a  check  on  the  accuracy  of  the  numerical 
scheme.  For  values  of  a  presented  in  Fig.  10  the  error  in  the  value  of  A  was  less  than  1%.  For 
Oi/Oo  oo,  p.  — >  1  and  equation  CB2)  becomes 

VT+VT 


-  f  g(t)  log  dt  =  g(s). 

Jo  V  s-t 


'0  V|  s-t , 

This  is  equivalent  to  the  linear  spring  problem  analyzed  fully  by  BUDIANSKY,  AMAZIGO,  and 
Evans  (198S)  and  ROSE  (1987). 


